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Abstract 

We analyze the eigenvalue density for the Laguerre and Jacobi /3-ensembles in the cases 
that the corresponding exponents are extensive. In particular, we obtain the asymptotic 
expansion up to terms o(l), in the large deviation regime outside the limiting interval 
of support. As found in recent studies of the large deviation density for the Gaussian 
/3-ensemble, and Laguerre /3-ensemble with fixed exponent, there is a scaling from this 
asymptotic expansion to the right tail asymptotics for the distribution of the largest 
eigenvalue at the soft edge. 

1 Introduction 

Let X be a p x N (p > N) random matrix with standard real Gaussian entries. It is well known 
in random matrix theory [26j H3] that the eigenvalue probability density function (PDF) of the 
corresponding covariance matrix W = X T X is proportional to 

N 

1=1 l<j<k<N 

with a = p — N + 1 — 2/0 and — 1. For general (11. ip is referred to as the Laguerre 
/3-ensemble and is to be denoted ME /3iA r(A /3a//2 e~ /3A/ ' 2 ). 

In multivariate statistics, thinking of X as a data matrix with successive columns storing 
the values of p measurements for each of the N observables, knowledge of (ll.ip allows for an 
analysis of questions such as "are the measured quantities in X correlated?" [[26]. In physics, the 
Laguerre /3-ensemble (II .ip with exponent N/2 appears as the eigenvalue PDF of the Wigner- 
Smith delay time matrix in mesoscopic transport problems [7] (see also [T3"t §3.3.1]). The 
parameter can then take on one of the three values = 1, 2 or 4 depending on the underlying 
time reversal symmetry. We remark too that the case — 2 of (11.11) corresponds to the 
eigenvalue PDF for the so-called Penner model in the theory of matrix models (see e.g. [28]). 
It also appears in the analysis of information capacity in the setting of wireless communication 
systems [30J. 

A peculiar feature of the Wigner-Smith time delay matrix example is that the parameter a 
in (11.11) is extensive, being proportional to the number of eigenvalues N. This is in fact what is 
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required in the analysis of the Penner model, and this same limit is of interest in the study of 
wireless communication systems; see references as cited above. One direction of study relating 
to fll.ip with a extensive is the large N asymptotic expansion of the moments 



N 



(1.2) 



1=1 



ME /J , JV (A' 3a / 2 e-/ 3A / 2 ) 



Recent references on this topic include [22], [231 [25]. These moments probe the non-oscillatory 
portion of the large N asymptotic expansion of the eigenvalue density in the domain of its 
leading support. Our interest in this paper is also in the large N expansion of the eigenvalue 
density with a extensive, but in the domain outside rather than inside its leading support. 

Physical motivation to pursue this problem comes from a very recent [16] experimental set 
up relating to fll.ip in the case (3 = 2, which in turn comes about by fll.ip corresponding to 
the eigenvalue PDF for matrices X'X, where X is a p x iV complex Gaussian matrix. In fact, 
the theoretical discussion in PS] tells us that the so-called round trip propagation matrix M 
determining the evolution of the complex electric field for a p x iV coupled array of high gain 
lasers can be decomposed in the form X'X for X a complex p x N matrix with approximately 
Gaussian entries. Furthermore, the output power is proportional to the largest eigenvalue of 
M. With hundreds of thousands of data points accessible in this way, it is possible to probe 
the large deviation portions of the largest eigenvalue PDF from these measurements, something 
that was done explicitly for the square case p = N(= 5). For the right tail the large deviation 
form of the largest eigenvalue PDF is, to leading order, equal to the density. This happens in 
the domain outside the leading order support, where the density takes on exponentially small 
values. 

There is another random matrix ensemble that allows for extensive exponents. This is the 
Jacobi /3-ensemble, specified by an eigenvalue PDF proportional to 



and referred to as ME ft j V (A /?ai/2 (l-A te/2 )). In the special case a 2 = 0, a x = p—N+1-2/0 and 
(3 = 1, 2 and 4 fll.3p gives the PDF of the non-zero squared singular values of an N x p sub-block 
of a (N + p) x (N + p) matrix from the circular orthogonal ensemble (/3 = 1), circular unitary 
ensemble ((3 = 2) and circular symplectic ensemble ((3 = 4) (see [321 §3.8.2]). The propagation 
of electron fluxes in a disordered mesoscopic wire is well described by scattering matrices from 
the circular ensembles. And an N x p sub-block then corresponds to the transmission matrix 
relating a state of p channels entering the left of the wire and a state of N channels exiting the 
right. The corresponding singular values then determine the conductance (see e.g. [5]). We 
will compute the asymptotic form of the spectral density for (11.31) with a\ and a 2 extensive, 
in the asymptotically small region outside the leading support. With (11.31) in the case (3 = 1 
occurring in mathematical statistics, and the largest eigenvalue having special significance in 
various statistical tests [19] . the fact that our results imply the right tail asymptotic form of the 
distribution of the largest eigenvalue thus provides us with motivation in an applied setting. 

Our task then is to compute the eigenvalue density for both the Laguerre and Jacobi (3- 
ensembles with extensive exponents, in the exponentially small regime outside the leading 



N 



nA? ai/2 (l-A,r 2/2 ][ |A*-A i 



< A; < 1 



(1.3) 
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support. In addition to the interest from the viewpoint of the problems discussed above, this 
complements the recent computation of the large deviation density for the Gaussian /3-ensemble 
[IH [5], and the Laguerre /3-ensemble with a fixed |14] . Here we will use the method of [14J 
based on writing the density as an average of the /3 moment of the characteristic polynomial, 
and then using a Gaussian fluctuation theorem to deduce the asymptotics. This is distinct from 
the method of [6J which is based on the loop equations for correlators. In Section 2 we review 
results on this problem already documented in the literature, and then present our general 
strategy. The details of our calculations, leading to the sought large deviation formulas, are 
given in Section 3. In the case of the Laguerre ensemble with weight \l 3aN / 2 +i 3 / 2 - 1 e -i 3X / 2 ^ where 
u L is given by (12. 9ft below and x G (0, (y/a + 1 — l) 2 ) U {{sfoc + 1 + l) 2 , oo), we obtain 



Np\ x) , N {Nx) r exp 



N0 



u L + 2 lo£ 



U + X 



a 



x exp { (1 - 3/3/2) log \u l 



x 



N{l + a)V 2 ( 2 



2tt 



(3N 



P/2 



2(1 + a) 1 / 2 
-(1-/3/2) log 

r(i + /3/2) 



2a log 



u 



x — a 



2x 1 / 2 {l + a) 1 / 2 
x{u h + x — 2 - 




:i-4) 



while for the Jacobi ensemble with weight A /3aiAr//2+/3 / 2 1 (1 — \y a2N / 2 +P/ 2 1 ) where u J is given 
by (I3.24p below and x G (0, ci) U (02, 1) with ci, C2 given by (12.31) . (12.41) . our final result reads 



P(i),.y(E)^^exp{iV£( log 



x- (Ci + C 2 )/2 + M J 



(ca-cO/2 



Oil log 



x — u 



a 2 log 



| V / (l-ci)(l-c 2 )-x + l + 



; v ^r + v ^)xV2 

x exp <( (1 - 3/3/2) log |w J | - (1 - /3/2) lo. 

xfciiri)' /2 r ( i + /3/2). 




(v/r^ + v/T^Xi-x)^ 

x(l-2;)(u J + 3;-(ci + C2)/2) 



2 + ai + a 2 



Ait 



(3N 



;i.s) 



The notation ~ in (ll.4p and (I1.5P means that LHS/RHS — > 1 in the appropriate limit. Both 
these asymptotic formulas are shown to exhibit a scaling to the asymptotic form of the right 
tail for the distribution of the largest eigenvalue at the soft edge. In the Appendix we critique a 
recent [20] claimed large deviation-type asymptotic expansion for the smallest eigenvalue PDF 
in the Laguerre /3-ensemble in the case that a is fixed. 
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2 Previous results and strategy 



2.1 Previous results 

For the Laguerre /3-ensemble with a = aN + 0(1) the leading support of the eigenvalue density, 
upon scaling the eigenvalues \i h-> NXi, is the interval (af , a|) with 

ai = V« + 1 - 1, a 2 = Va + 1 + 1. (2.1) 

Moreover, the corresponding leading eigenvalue density is given by the so-called Marcenko- 
Pastur law 

N 



Np^(Nx) ~ y(a| - x) (x - af), a{<x<a 2 2 (2.2) 

(see e.g. [131 ec L- (3.61) with x 2 h» a;]). Thus for a > the smallest eigenvalue is bounded away 
from the origin with a density profile that vanishes like a square root (a soft edge), while for 
a = the smallest eigenvalue is part of an accumulation of eigenvalues about the origin with 
a density profile diverging like a reciprocal of a square root (a hard edge). 

A similar effect is exhibited by the Jacobi /3-ensemble f ll.3p with a\ = a\N + 0(1), a 2 = 
a 2 N + 0(1), except now there is no need to scale the eigenvalues. Thus the leading eigenvalue 
support is (ci,c 2 ) C [0,1] with 



2{a\ - g\ 
(ai + a 2 + 2f 



2(0, +<*-!)= ,: m :" 2 ' (2.3) 



(2c 1 -l)(2c 2 -l) = r M±j) _ L (2 . 4) 
(ax + a 2 + 2) 2 

Notice that (12. 3ft and ( 12. 4ft are unchanged by 

(ci, c 2 , ai, a 2 ) (1 - ci, 1 - c 2 , a 2 , a x ). (2.5) 

The corresponding eigenvalue density is given by 

2 + gi + a 2 y/(x - ci)(c 2 - x) 
Pd)^)-^^^ ^x) (2 ' 6) 

(see e.g. |T3j eq. (3.77) with y t-^ 2y — 1]). 

We are interested in the eigenvalue density outside of [af,a|] for the Laguerre /3-ensemble, 
and outside of [c\, c 2 ] for the Jacobi /3-ensemble. In the former case there are some prior results 
from the existing literature, which we will now summarize. But before doing so, we make note 
of the work [23]. This relates to the case a = in the Laguerre ensemble, and so strictly 
speaking is not relevant to our setting of extensive exponents. However, the leading term of the 
asymptotic expansion of the density does not contain any information on the exponent a for a 
of order unity, and so is expected to correspond to the case a = 0. Actually the concern of [23J 
is not the density in the exponentially small region as such, but the large deviation asymptotic 
form of the distribution of the largest eigenvalue P^b(s). However, it is a simple result that for 
x in this regime, the two are identical [131 eq. (14.136)], 

jkjiNx) ~ fa{Nx). (2.7) 
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The leading large N form of (12.71) for x > 4 in the case a = as calculated in [23J is given by 
setting a = in (12 .8p below. 

There are at least three existing works containing explicit asymptotic formulas for (Nx) 
(or equivalently, according to (12.71) . p^^Nx)) outside of [af,a|] for general a. In order of 
publication date the first [US, ec L- (4-83)] is the leading large N form 



exp 



(3N 



-u + a Ioe 



a(a + u L — x) — 2x 



2(a + iy/ 2 x 



(2 + a) log 



u + x 



a 



2(a + l)V2 



where 



u 



l ._ 



yj (x — a\)(x — a 2 ,), x e (a^oo) 
(x — af)(x — a|), xG(0,af). 



+ o(iV) 
(2.8) 

(2.9) 



The second [20] is presented as the leading asymptotic form for p^^(Nx) in the region x G 
[0, a 2 ]. In terms of our notation the result of [20J reads 



pL } (iVx) = exp (-/3iV$ min ((1 - v^TT) 2 - x) + o(JV)) 



(2.10) 



where 



-Wx x + 4Va+T - - log 1 



a 



(i - v^+iy 



+ 2 log 



\/x + + 1 — \fx 



+ a log 1 + 



\/Va + 1 

2^ / a/x + 4a/o + 1 - v/i' 



VaT+T- 1 



4 v / aTT 



(2.11) 



Elementary algebraic manipulation involving the simple identities 



a log 



a(a + u — x) — 2x 



2(a + l) 1 /2 



a log 



u L + x-2- a 



2(a+ iy/ 2 x 



2a log 



a + u + x 



2y/a + ly/x 



—2a log 



a — u + x 



2y/a + ly/x 



(2.12) 



shows that (12. 8 p and (I2.10p are in fact identical expressions. 

The third and most recent result [T] was derived for the case (3 = 2 only, but in addition 
to the leading exponential term also gives the algebraic correction. Thus, after making the 
correspondences A = a,(, = a\/(a + l),£i = a 2 ,/ '(a + l),x h> x/(a + 1) we read off from [TJ 
that for exponent a = aN 



pU Nx ) 



(a 



+ l) 1 ' 2 



2irN (x - a\)(x - al) 



-Nk(x) 



(2.13) 
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where 



k(x) = —u L + a\a 2 arcosh 
Recalling that for x > 1 



I - I I 1 -u 1 

x 2 I 5? I" ^| 
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i) 



a 2 + a 2 . 



arcosh 



— X 



(2.14) 



arcoshx = log (a; + V x 2 — 1) 
shows that the leading (exponential) term of (12.131) coincides with the (3 = 2 case of (12. 8p . 

2.2 Strategy 



We seek to extend ( 12. 8ft to the same order as the (3 = 2 result (12.131) . and thus give the 
asymptotic form of pju(Nx) for x G [0, af ) U (a|, oo) up to terms o(l) as iV — )■ oo. We seek 
too an analogous expansion in the Jacobi case. Our method relies on the fact that the density 
can be written in terms of the /3-moment of the corresponding characteristic polynomial, and 
as such is an application of the strategy used recently in the case of the Gaussian /3-ensemble, 
and the Laguerre /3-ensemble with a fixed |14j . 

We have referred to the PDFs JTTJ and by ME p>N (w) with w(X) = X a ^ 2 e~^ x / 2 and 
w(X) = A ai/3 / 2 (l — A)" 2 ^ 2 respectively. Let us denote the corresponding normalizations by 
C^tvM- Then it follows from the definitions (see e.g. [TJJ eq. (2.11)]) that 



(iv + i^M 

P(i),v+i(^) = — ~ r-i — 



A/ 



(2.15) 



MEu,jv(i«) 



The key observation is that the average in (I2.15P can be interpreted as a particular case of the 
characteristic function for the linear statistic V(x) = YliLi 1°§ \ x ~ But for a large class of 
weights w the latter is expected to be a Gaussian (see e.g. [131 §14.4]), and this implies, with 
x outside the interval of support, 



N 



~ e^M e (M»)) 2 /2 



(2.16) 



ii=i 



ME PtN (w) 

where, with v(t) := log \x — t\ and (d\, d 2 ) the limiting interval of support 

\ P(i),N(t)v(t)dt 
'di 



(2.17) 



dt\- 



d - v'{t 2 )({d 2 -t 2 ){t 2 -d x )) 
dt 2 1 



1/2 



2 r f d x + d 2 d 2 -d x \ 

a,k = — I v\ 1 cos0 J coskv da. (2.18) 

Wo I 2 2 J 
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Theorems relating to the asymptotic form of the characteristic function for linear statistics 
in random matrix theory were first proved in [17J. However, the technical assumptions therein 
did not permit the linear statistic corresponding to v(t) = \og\x — t\ due to its singularity 
at x — t. But for x outside the eigenvalue support (I2.18P is well defined and on physical 
grounds it is expected that (" 12 . 16[) will still apply [131 §14-4]. Fortunately, very recently [5], 
under the conditions that the eigenvalue support is a single interval, and with both endpoints 
soft edges (or a hard edge with parameter a = 0) the validity of (I2.16P for any v analytic in the 
neighbourhood of [di,<i 2 ] has been rigorously established. In our setting of the Laguerre and 
Jacobi ensembles with extensive parameters both end points of the eigenvalue support are soft 
edges, so (" 12 . 161) thus follows from [5]. 

Our immediate task then is to compute the explicit form of (I2.17P and ( 12 . 1 8j) . We will 
perform this task in the next section. Substituting (I2.16P in ( I2.15P then gives the sought 
asymptotic expansions of the densities. 

At a technical level, with the parameters being extensive, and the occurrence of both N + 1 
and N in ( I2.15p . there is some advantage in first manipulating ( I2.15P before applying ( I2.16p . 
Thus in the Laguerre case the weight required in (I2.15P is 

w f x \ _ x a(N+l)/3/2+l3/2-l e ~{N+l)/3x/2_ ^.19) 

With this w we observe 

where w(x) = x aN P/ 2+ P/ 2 ~ 1 e~ Nl3x / 2 , and we observe too that 

C/3,n[w] _ 1 



(2.21) 



Similarly, in the Jacobi case, the weight required in (I2.15P is 

w f x j = x ax(N+l)p/2+p/2-lQ _ x ^a 2 (N+l)p/2+p/2-l (2.22) 

and we have 

f[\x-\A =^^(flX^ / \l-X l r^\x-\A , (2.23) 

where w(x) = x ^P/2+p/2-i^ _ x }a 2 Np/2+p/2-i^ together with 

Cp, N [w] _ 1 



(2.24) 



Our immediate task then is to compute the explicit form of ( 12.1 7p and (I2.18P for the averages 
on the RHS's of (12T20|) and (12~23|) with 

„(t\ = I log \ X ~ *l + f lo S* ~ I*' Laguerre 
K ) \ \og\x-t\ + aiogt+^l g|l-t|, Jacobi, y > 
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and with 



v(f) _ J f logt - Laguerre 

UW ~ \ f logt + f log |1 — t|, Jacobi. [Z ' M) 

For future purposes, especially in labellings, it will be convenient to refer to ( 12.251) and ( 12.261) 
as choices (1) and (2) respectively. We will proceed with required computations in the next 
section. 



3 The asymptotic forms 

3.1 The Laguerre ^-ensemble 

For the weight w in (I2.20p it is known from [TT] that 



N I 

P(i),N(t) = — \](al -t){t- a\) 



2vrt 



+ 




^(t-a?) + ^(t-al) 



V(t-a?)(al-t) 



+ 



71 



N 



(3.1) 



We remark that the precise form of the 0(1) term in (13.11) relies crucially on the 0(1) portion 
of the exponent in w equalling (3/2 — 1. Note that the latter is precisely the 0(1) portion of 
the exponent appearing in (II. ip in relation to covariance matrices. We remark too that the 
expansion ( 13. ip ignores possible oscillatory terms which appear at next order [9], as these do 
not contribute to the next correction in (I2.17p . 

As is consistent with the labelling noted below (I2.26p . let ( I2.17P in the case p(i),v(t) is 
given by ( 13. ip . and v (t) by the Laguerre cases of ( I2.25P and ( I2.26p . be denoted /i^ and 
respectively. Then fJ^^ — l^f 2 ^ is equal to (I2.17P with pn\sit) given by ( 13. ip and 
v (t) = log | a; — t\. It can be evaluated by making use of the following integral evaluations. 

Proposition 3.1 Let u h be given by \2.9\) . For x G (0, a\) U (arj, oo) we have 



2tn 



(t — af)(a| — t) log \x — t\ dt 



x 



a 



u L - 2 + a log 



a (a + u 



x) 



2x 



and 



log \x — 1\ 



y/(a%-t)(t-a*) 



2x 2 



dt 



+ (2 + a) log 



u L + X 



a: 



(3.2) 



1 



a 2 - a 2 . T 

log — + log - x ± vx? — 1 

6 2 2 



log 



u h + X 



a 



(3.3) 



where in the first equality of 113. 3\) we take + (— ) according to x G (a 2 ,, oo) (x G (0, a 2 )) and 

2 



x 



x — 



a\ + a 2 



(3.4) 
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Proof: The result ( 13 .2p is derived in [131 Exercises 14.4 q.6(i)] for Re(x) > a 2 ,; its value for 
x G (0, a 2 ) follows by analytic continuation. The first equality in (I3.3P is a simple corollary of 
the integral evaluation 



1 f 1 log \x — t\ 



1 



dt = log ( - (x + V x 2 — 1) 



x > 1, 



the derivation of which can be found, for example, in [HI eq. (3.3)]. The second equality in 
(13. 3p now follows by noting 



Qir% Qj-\ 



Vx 2 -1= u L . 



(3.5) 



□ 

Proposition 13.11 implies that (I2.17P in the case p(i),jv(t) is given by (13. ip and v(t) by the 
Laguerre case of ( I2.25P has the large N evaluation 



L,(l) _ L,(2) 



N 



+ (2 + a) log 



x — a — u L — 2 + a log 

u L + x — 2 — a 



a(a + u — x) — 2x 



1 1 

^"2 



2x 2 
log \u L \ — lo. 



u + x — 2 — a 



+ 



1 

N 
(3.6) 



It remains to evaluate (12.181) with v{t) given by the Laguerre cases of (12.251) and (I2.26p . 
Consider the former case first. Making use of the second equality in ( I2.18P shows we must 
evaluate 



n _ n W , _(2) , (3) 



k = l,2, 



where 



with x given by (13 .411 . and 



7T 



log 1 - 



cos 9 



X 



cos fc^ dd, 



(2) « (1) 

«1 =2 a k 



x=0 



(3) 
«fc = 

From [HI Lemma 2] we know that 

cos 9 s 



at — a 



2 />7T 



2tt 



cos 9 cos ci^. 



log 1- 



X 



log(l - ^e ie )(l - v x e- 10 ) - log(l + v 2 x ) 



-ie\ 



where v x is such that \v x \ < 1 and has the explicit value 



x-(x 2 - l) 1 / 2 , 
£ + (x 2 -1) 1/2 , a;e(0,a 2 ). 



(3.7) 
(3.8) 

(3.9) 
(3.10) 

(3.11) 

(3.12) 
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We see immediately from ( 13 .lip that 



k 

According to CT) . we see from this, (1512]) . (ET3j) and (EE} that 

,( 2 ) - 1 



«?> = - 2 4- (3.13) 



K y/a + 1 



And in relation to (13 . 10 j) . an elementary calculation shows 



af = { " V "^ 1 ' Z =L - (3-15) 



--v/cH-T, fc = 1, 
0, otherwise. 

Substituting these explicit forms in ( 13 .7p we see 

^ ka\ =a + 1 + 2Va + 1 (2z/ x + cw/ ) + ^ - (2z/ + az/£) 2 
fc=i fe=i 

= - a + 1 + 2(x - (a + 2) - u L ) - 41og(l - ^) 

- 4alog(l - z/^z/o) - « 2 log(l - vl). (3.16) 

With t>(t) given by the Laguerre case of ( I2.26H . in the notation of ( 13. 1\ we have = 
+ af\ The first line of (13. 161) then remains valid but with v x = 0, so in this case we have 

oo 

^ka\ = -a + 1 -a 2 \og{l-vl). (3.17) 

k=l 

Recalling now ("12. 18[) and using a superscript notation analogous to that used in (I3.6P we 
read off from <KT^ and (IBTTTj) that 



(<7 2}L,(i) _ f<7 2^L,( 2 ) = }_(^ 2{x _ ( a + 2 ) - M L ) - 41og(l - ^) - 4alog(l - ^ )). (3.18) 



Use of (13.121) and (13. 14j) allows us to compute that 

\x-(a + 2) + u L \ 



log(l - v x ) = \og\u I - log 
log(l - u x uo) = log (3-19) 

Substituting in (13.181) then gives the explicit evaluation of the variance 

^ , o^ ..m 2 i^,„.l, , 2 la;-(a + 2) + M L | 



((T 2 )L ,(1) _ (^,(2) = _ (x _ (Q + 2) _ U L) _ lQg | U L| + lQ? 



2a \x + a-u L \ 
/3 2(a + l) 
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According to ( 12 . 16[) the asymptotic form of the product of the averages in (I2.20p and (I2.2ip 
is determined, up to terms vanishing as iV — >■ oo, by (13. 6p and (I3.18p . allowing us to read off 
from these that 



n \ x - Xi 







J=l 



ME^ i ^(A«( iv + 1 )' 3 / 2 +/ 3 / 2 - 1 e-( JV + 1 )' 3A / 2 ) 



(N + l)8 i , 
~ exp < ( x - u L - a - 2 + 2 log 

TV— »00 I 2 



u + x — 2 — a 



2a log 



u — x — a 



2(1 + a 



x exp <^ (1 - 3/3/2) log - (1 - 0/2) log 



U + X 



a 



(3.21) 



Although in the derivation x has been assumed fixed, taking x — > oo shows that the RHS tends 
to x ° , which is the leading form of the LHS. This suggests that in fact (13.211) is valid uniformly 
for x G (0-2,00). We remark that the average of YI1L1 \ x ~ f° r the Gaussian /3-ensemble, 
and the Laguerre /3-ensemble with exponent a fixed, was shown to exhibit the same property 



The announced formula (11.41) follows by substituting (13.211) in (12.151) together with the 
asymptotic formula for the ratio of normalizations 

2 V /2 



C/3,N [A 



a(iV+l)/3/2+/3/2-l -(Af+l)/3A/2 



1 



C/3,7V+1 [\ a (N+^/2+P/2-l e -(N+l)0X/2+fi/2-lj 



J{N+l)(l+a/2) 



2tt \(3N J 
x r(l + /3/2)(l + a) 



-(jV+l)(l+a)/3/2+l/2 



(a consequence of the explicit gamma function evaluation [T3J Prop. 4.7.3] and Stirling's for- 
mula). This formula is in precise agreement with the result (12.131) from pQ. With a = and 
x scaled x \-> 4x (II. 4ft reduces to the a — 1/2 — 1/0 case of the asymptotic formula for the 
density p^(4iVx), x > 1, in the Laguerre ensemble with weight w(\) = \l 3a / 2 e~ l3X / 2 given re- 



cently in [H] (the significance of the choice a = 1/2 — 1/0 is that the Laguerre ensemble then 
reduces to the chiral ensemble with weight w(\) = e~^ A / 2 ; see e.g. the discussion in [15j §3.2]). 
For = 1, 2 or 4 the validity of (11.4j) can readily be illustrated numerically by comparing its 
evaluation against that of the exact expressions for the density known in terms of Laguerre 
polynomials (see e.g. [131 Ch. 5& 6]). For definiteness consider (3 = 2, for which 



and (3 



T(N+1) 
T(N + a) 

1, N even, for which 



x e 



x 



- L a N _^ 



x 



d 

dx 



L jy (X , 



u (A)=A Q / 2 - 1 e- A / 2 



0=2 



_ x (a-l)/2 e -x/2 L a 



N-l 



X) 



X- 



r(A0 



sgn(x 



t)L a N _ 2 (t)tW 2 e-*/ 2 dt. 



4r(a + N - 1) J(l 

In these expressions set a = <x/V, and denote the corresponding asymptotic form of as implied 

L, asym / 



by (H3D by p^ asym (Nx). We now perform a numerical computation of the ratio 

r\ N {x-o) 



P^ Sym {Nx) 
p h {l) {Nx) 
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0.1 


6 


10 


0.1 


6 


10 


6 


1.670 


4.318 


1.122 


1.072 


1.205 


1.083 


12 


1.357 


2.978 


1.062 


1.049 


1.101 


1.041 


18 


1.246 


2.460 


1.041 


1.039 


1.065 


1.027 


24 


1.189 


1.175 


1.031 


1.027 


1.048 


1.020 


30 


1.153 


1.991 


1.025 


1.023 


1.039 


1.016 



Table 1: The ratios r\ N (x; 1) (first three inner columns) and r\ N (x\ 1) (final three columns) 
for varying iV (rows) and given x (columns). The eigenvalue support is (0.171,5.828) 



The obtained results for given x and varying N are given in Table [TJ 



3.2 The Jacobi /^-ensemble 

The derivation of the first non-oscillatory correction term for the Gaussian /3-ensemble density 
given in [131 §14.3] implies that in the case of the Jacobi /3-ensemble with weight as specified 
below (|1.5p . (12. 6p can be extended to read 



p {1)iN {t) = N 



2 + ai + a 2 V(*-Ci)(c2 -t) 



2tt 



t(l-t) 



+ 



1 1 

^"2 



1 



b(t - Cl ) + h(t - c 2 ) - \— = — 

2 2 * \f{t-c x ){c 2 -t) 



O 



N 



(3.22) 



(cf. the first asymptotic equality in (13. ip ). In a notation analogous to that used in ( 13.61) . to 
compute fi J ^ — requires computing ( 12 . 1 TP with p(i)^(t) given by ( 13.221) and v(t) = 

log \x — 1\. 

To carry out this task, let the first term in (I3.22j) be denoted by JV^, oo( x )- We know from 
[31] (see also [27], [25]) that 



where 



(-2 



Cl 



P(i),oo\ x ) , ot\ a 2 

— dx = h 

w — x w 1 — w 



- (oti + a + 2) 



u 



w(l — w) 



(3.23) 



it' 



a/( w - Cj)(w - C 2 ), > c 2 

-a/ (w - ci)(w - C 2 ), IW < Ci. 



(3.24) 



The anti-derivative of the final term can be computed using computer algebra. Choosing the 
integration constant so that as w — > oo the constant term vanishes tells us that 



f° 2 1 ( 

2 J P(i),oo( x ) log \x -X\dX = + a 2 + 2) log - ^ 



Cl + C 2 j 
h li 



+«i log 



i?(x;ci,c 2 ) 



x- 



+ a 2 log 



R(l — x; 1 — ci, 1 — c 2 ) 



;i-x) 2 (^r^7- v/r^) 



(3.25) 



12 



where 

R(x; ci, c 2 ) := x(c\ + c 2 ) - 2cic 2 - 2 A /cic 2 u J . 
In fact simple manipulation shows the RHS of (I3.25p can alternatively be written 



(3.26) 



2 log 



C l + C 2 , J 
- X h U 

2 V 2 



— 2a i log 



2|v/CiC 2 + x - u J 

(v^r+v^) 2 



- 2a 2 log 



2|V(1- Ci )(1-c 2 ) + 1-x + m j | 



(3.27) 



The integral required by the second term in (I3.22p to compute (12.171) with v(t) = log \x — 1\ 
is given by (13. 3p . Thus we have 



J,(l) _ J, (2) 



N I** p\ 1)t00 (X)\og\x-X\dX 

Jci 

(~0 ~ 2) (2 ^ I ( X ~ Cl ) - C 2) I - log ^ 



C l + C 2 , J 
X h W 



(3.28) 



We now turn our attention to the calculation of the variance. This task is very similar to 
that already detailed in the Laguerre case, and so we will be content with presenting the final 
result only. Thus we find 



2U(1) 



2\J,(2) 



2 . j. 2 

— log \u I + — k>£ 



8 



B 



C l + C 2 , J 

21 1 2— +U 



2ai . 2\JciC 2 + x - u J 
— log- 



2a 2] 2|v/(l- Cl )(l-c 2 )-(x-l-n J ) 
log — 



B 



S2") 2 

(3.29) 



Substituting ( 13. 28ft and ( 13. 29ft in ( 12. 16ft and recalling (I2.20p allows us to deduce the asymp- 
totic expansion 



N 



n i x - a ' 



~ exp <(N + 1)8 log 

TV— »oo 



1 



Cl + c 2 j 
X h IT 



-ai log 



2|v/cic 2 + x - u 3 



«2 log 



2| v /(l-c 1 )(l-c 2 )-x + l + n J 
( v ^7^ +v ^^)2 

1 



xexp < (1 — 3/3/2) log |w J | - (1 - 8/2) log - \u J + x - (c x + c 2 )/2 



(3.30) 
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0.8 


0.85 


0.9 


0.8 


0.85 


0.9 


6 


1.866 


1.224 


1.115 


1.045 


1.055 


1.055 


12 


1.487 


1.116 


1.058 


1.030 


1.040 


1.028 


18 


1.345 


1.079 


1.038 


1.033 


1.028 


1.018 


24 


1.269 


1.059 


1.029 


1.035 


1.021 


1.014 


30 


1.221 


1.048 


1.023 


1.035 


1.017 


1.011 



Table 2: The ratios r\ N (x; 5, 5) (first three inner columns) and r\ N (x; 5, 5) (final three columns) 
for varying N (rows) and given x (columns). The eigenvalue support is (0.223,0.776) 



From this, (I2.15P and the asymptotic formula for the ratio of normalizations 

Q ^ax{N+i)/2+l3/2-l(i _ ^0a 2 (JV+l)/2+0/2-l) 
C^ iA r+i(A' 3a i( JV+1 )/ 2 +/ 3 / 2 ^ 1 (l - A)0 a 2(W+l)/2+0/2-l) 

T(l + /3/2) / 2 \^/2 ((l + a 1 )(l + a 2 )(l + a 1 + a 2 )) 1 / 2 
2tt \(3n) (2 + ai + a 2 yi 2+l 

111 

x exp |/?(iV + 1) ^2 log(ai + a 2 + 2) - - bg(a a + 1) - - log(a 2 + 1) - - log(ai + a 2 + 1)) 
+f3(N + l) ai ( log(ax + a 2 + 2) - ~ log(ai + 1) - ~ log(ai + a 2 + 1)) 
+P(N + l)a 2 ( log(ai + ol 2 + 2) - - log(a 2 + 1) - - log(ai + a 2 + 1)) } 

(deduced from the Selberg integral [T3l Ch. 4] and Stirling's formula), and the inter-relations 

2(l + a 1 ) 1 / 2 (l + a 1 + « 2 ) 1 /2 



\A - C 2 + y/1 - Ci 



C 2 - Ci 



(2 + «i + a 2 ) 
2(l + a 2 ) 1 / 2 (l + a 1 + Q2 ) 1 / 2 

(2 + «i + a 2 ) 
4(1 + aijV^l + a 2 ) 1 / 2 (l + ai + a 2 )V 2 
(2 + a! + a 2 ) 2 



(deduced from (Q and (Q) we obtain (Oil . 

As for fll.4p . exact formulas for the density at /9 = 1, 2 and 4 (see e.g. [131 Ch. 5&6]), now 
involving Jacobi polynomials, can be used to check ( 11. 51) . Explicitly, consider the Jacobi weight 
specified above ( II. 5p . and let Pny^ix) denote the asymptotic density for x G (0,ci) U (c 2 , 1). 
The ratio 

can then readily be computed numerically for 01,2 and 4. As seen in Table [2j results similar 
to those presented in Table [1] are obtained. 
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3.3 Soft edge scaling 

For the Laguerre /3-ensemble with a = aN there is a scaling about the smallest and largest 
eigenvalues such that the neighbouring eigenvalues have a spacing 0(1). Thus about the largest 
eigenvalue this occurs in the variable X when [18] 



NX = N(l + Va + l) 2 + N 1/3 (Va + l + 1) ( 1 + . 1 ) ^ X (3.31) 

V \/ni. 4- 1 / 



and for the smallest eigenvalue when [2] 



NX = N(l - V^Tlf - N 1/3 (V^TT - 1) ( 1 - . 1 ) X. (3.32) 

v y/a + 1/ 

We note that with t± := (vaT+1 ± l) 2 , these can be written 

^^«" !,3 ((dS^' (3 - 33) 

where the +(— ) sign is to be taken about the largest (smallest) eigenvalue. Moreover, the 
distributions of the largest and smallest eigenvalues have well defined limits in terms of X, 
independent of a. These are the soft edge distributions, common — according to the principle 
of universality [21] - to a large class of random matrix ensembles at the spectrum edge, and 
in particular the Gaussian /3-ensemble and the (soft edge of) the Laguerre /3-ensemble with 
parameter a fixed. For these latter ensembles it is known [101 HH E] that the distribution of the 
soft edge scaled largest eigenvalue, p soit (X) say, has the large X form 

n^(X) i r(l + /V2)e^ 3/2 /3 

Replacing x in (jl.4p by the RHS of (I3.33p . then taking iV — > oo we reclaim f)3.34p . 

There is an analogous scaling in the case of the Jacobi /3-ensemble with a\ = a%N, a 2 = a 2 N. 
About the largest (+ sign) and smallest (— sign) eigenvalues this reads [191 S] 

\ _ +J -U AT- 2/3 i^ii 1 ~ 4)) 2/S /o op-X 

A " t±±iV (r(9(l-r)(l-e))V6' ^ 



where t\ — C2, tL — c\ and 



9 := -, r : = 



cki + a 2 + 2 q;i + a 2 + 2 

Replacing x in (jl.5p by the RHS of (13.351) and taking — > oo again reclaims (I3.34p 
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A Appendix 



Consider the Laguerre /3-ensemble with weight A a/3//2 e~^ A//2 and suppose the exponent a is fixed. 
The eigenvalue support then extends to the origin, so there is now no region to the left of 
the interval of support for which the density is exponentially small. As first noted in [8], the 
appropriate large deviation quantity is now Ep^(0; (0, s)) — the probability that the interval 
(0, s) is free of eigenvalues. Computation of the large deviation asymptotics (using the Dyson 
log-gas method), setting s h> X/4N, and taking iV — > oo reproduces the X — > oo asymptotics 
of the hard edge scaled gap probability 

4 ard (0;(0,X)):= lim E P)N (0, (0, X/AN)). 

^ N— s-oo 

According to [12] (see also the recent works [131 ec L- (13.52)] and [29J) the latter reads 



lo gB r(0;(0,A-)) x ^-/3(f -a^ + (^ + JL) log A^ + (l)) (A.l) 

(for a G Z> the explicit form of the constant term is also known; see e.g. [131 below eq. (13.52)]). 

There has been a recent work relating to the asymptotic form (IA.1I) |20j from the large 
deviation viewpoint. There (eq. (21) of [20J with z G [a 2 /4, oo)) an expression functionally 
simpler than ( lA.lj) is presented. However when compared against (lA.lj) it is seen to differ in 



the coefficient of log A 1//2 , and differ too in the constant term in the cases the latter is known 
exactly (the large Nt region of the graph in Figure 2 of [20] also reveals this discrepancy) . The 
formula eq. (21) of [20] must therefore be regarded as an approximation rather than an exact 
result. 
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